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Abstract 

The mean-field dynamics of a collection of stochastic 
agents with local versus nonlocal interactions is studied 
Via analytically soluble models. The nonlocal interactions 
result from a barycentric modulation of the observation 
range of the agents. Our modeling framework is based 
,on a discrete two- velocity Boltzmann dynamics which can 
be analytically discussed. Depending on the span and the 
modulation of the interaction range, we analytically ob- 
serve a transition from a purely diffusive regime without 
definite pattern to a flocking evolution represented by a 
solitary wave traveling with constant velocity. 

Keywords: Self-organized systems, Transport processes, 
Kinetic theory. 

1 Introduction 

Since several decades, the fascination for flocking dynam- 
ics and specifically, the detailed mechanisms which induce 
a collection of interacting stochastic agents to exhibit an 
emergent collective behavior, stimulates a fruitful activity 
on both experimental and modeling sides. The recent pio- 
neering and truly engrossing contributions from T. Vicsck 
et al. HZ] and F. Cucker and S. Smale [TT] unveil, from 
a simple and synthetic modeling point of view, some of 
the basic features underlying the formation of robust col- 
lective motions. As intuitively expected, a central role is 
played by the interactions connecting the agents. Both the 
interaction range and the interaction strength are control 
parameters which, for critical values, may trigger a flocking 
(phase) transition. By a flocking pattern, we understand 
the self-organized capability of agents to create, via their 
mutual interactions, a persistent probability density tak- 
ing the form of a traveling wave. Hcuristically, with strong 
enough interactions, the noise-induced detuning tendency, 



can be overcome and ultimately produce an emergent syn- 
chronized evolution. In the models introduced by Vicsck 
et al. and by Cucker and Smale, it is indeed the interaction 
range that directly tunes the self-organization capability. 
Alternatively, recent observations for birds demonstrate 
that a fixed critical number of neighbors rather than a 
pure metric spatial range is used by birds in order to flock 
together to platoons [SJIE]- This shows that the interaction 
range itself may depend on the instantaneous distribution 
of agents. Here, we contribute to this widely open topic by 
proposing a simple class of models, where the interaction 
range depends in an effective way on the agent distribu- 
tion. 

Societies of agents can be composed either of dynamically 
homogeneous or heterogeneous individuals and each case 
requires a drastically different modeling approach. Re- 
cent analytical, while impressive, results have been de- 
rived for heterogeneous agents like ranked-based interact- 
ing Brownian motions [21 HJ [M] or bucket brigades dy- 
namics [TJ 3] . We here focus on homogeneous populations 
of agents which offer a wide potential for analytical ap- 
proaches. In an homogeneous population, any randomly 
selected individual is likely to be a dynamical representa- 
tive of any other fellow of the society. This very basic fea- 
ture together with the fact that for large populations, the 
relative importance of fluctuations diminish, enable, in the 
thermodynamic limit N — > oo, to adopt a purely macro- 
scopic description. This assumes that the normalized dis- 
tribution of agents can be "hydrodynamically" described 
by a probability density which solves a deterministic, yet 
nonlinear, evolution equation. We call the resulting equa- 
tion the mean-field (MF) equation. 

The class of MF models introduced in this contribution 
shows flocking transitions for large enough interaction 
ranges and belongs to the rare instances of exactly solv- 
able interacting agent dynamics. The price to be paid 
for such an analytical objective is to conveniently reduce 
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the individual agents' state space and decision space. As 
a minimal model, we shall consider random two velocity 
models involving agents traveling on the real line with two 
discrete velocities, say {v±} with w_ < v+. The effec- 
tive agent autonomous decisions will consist in selecting, 
at non-homogeneous Poisson random times, one out of 
the two possible velocities. The velocity updating process 
will, via density-dependent parameters of the underlying 
Poisson switching times, depend on the observation range. 
This naturally leads to nonlinear Boltzmann-like dynam- 
ics, which is commonly adopted for traffic flows model- 
ing HH [18]. The observation range, which will be the 
key control parameter for our discussion, may depend on 
the distribution of the whole agent society via a barycen- 
tric modulation. Such a barycentric modulation offers, in 
an effective manner, the possibility to model agent dynam- 
ics with configuration-dependent interaction ranges. This 
approach enables us to construct several solvable nonlin- 
ear models for which critical interaction modulations - re- 
quired for flocking - can be obtained explicitly. 
Although different from our present contribution, related 
MF models exist in the literature. Exactly solvable MF 
models have been proposed in different contexts including 
portfolio theory [TJ], coupled phase oscillators [7J, traf- 
fic dynamics [5] \W\, self-propelled organisms [3B], etc. 
Related MF models with barycentric self-interactions are 
studied [2] and [ID] . Moreover, the class of models dis- 
cussed in the present paper is amenable to a discrete ve- 
locity Boltzmann equation of the Ruijgrok-Wu type [25] 
and can be hence analytically discussed. As originally dis- 
cussed in [5T][52], f° r large populations of agents, our dis- 
crete velocities Boltzmann dynamics can be made to con- 
verge towards the generalized Burgers equation for which 
exact results and stability issues for traveling waves can 
be discussed explicitly, [2"3"] . 

2 Two velocity agent model 

Consider a collection A composed of N autonomous agents 
which are in a migration process on the real line M. At 
any time t S R + , we assume that the complete popu- 
lation is composed of two types of agents A.+ and A— , 
(A = A+ U A- ) , characterized by two associated migra- 
tion velocities v + and w_ on R. Agents traveling with ve- 
locity v + (resp. V-) belong to A+ (rcsp. A-). The time- 
dependent positions of the agents Xf.(t), k = 1,2, ...,N, 
can be written as a set of coupled stochastic differential 
equations (SDEs): 

X fc (t)=J Jk (X(t)), k = l,2,...,N, (1) 

with X(t) = (Xi(t), ...,X N (t)) and where 4(X(i)) stands 
for a two-states Markov process with state space £1 := 
{v+,V-}. The associated transition rates of agent k, 
a = a(Xk(t), X(i)) > from v+ to v_ and /3 = 
P{Xk(t), X(i)) > from to v + are here state- 
dependent. For a given configuration X(t) at time t, they 
represent the inverse of the average sojourn times of agent 



k in the velocity states. This dependency of the switching 
rates on the agent population effectively allows agent k to 
observe his/her environment and to react accordingly. In 
the sequel, we consider large populations of homogeneous 
agents, i.e. N — > oo, so that a MF description of the 
dynamics holds. Using a phenomenological MF represen- 
tation, we shall write P+{x,t \ xq) (resp. P-(x,t | xo)) for 
the conditional probability density to find agents at posi- 
tion x at time t with velocity v+ (rcsp. i>_), knowing that 
at time t = the density was given by the initial distri- 
bution p+{x) (resp. p-(x)). Due to the assumed Markov 
character of the transitions, the nonlinear evolution for 
P+{x,t | xo) and P—(x,t \ xq) can be written as a discrete 
Bolzmann-type equation (n.b. we drop the xo-dependence 
for typographic ease): 

P±(x,t)+v±d x P±(x,t) = 

T a{x, t)P+ {x, t) ± (3{x, t)P- {x, t), (2) 
with interaction kernels: 

rx+V 

a(x,t)=a- g[z- {X(t))]P-(z,t)dz, (3) 

Jx-U 
rX+V 

p(x,t)=p+ g[z-{X(t))]P + (z,t)dz, (4) 

Jx-U 

and where now a and /3 are two positive constants, g{x) > 
is a smooth function and integrable with respect to the 
P+ and P- probability densities and where 

(X(t))= [ x[P + {x,t)+P_{x,t)}dx (5) 

is the barycenter. According to eqs. ([2J5]) , agents modify 
their velocity by the following dynamic rules: 

(i) Agents with velocity w + . An agent at position x 
and with velocity v + changes spontaneously to V- 
with constant rate a. This agent is allowed to ob- 
serve the environment in the observation interval O = 
[x — U, x + V] . The a-rate is reduced in a weighted 
proportion to the number of v— agents present in O. 
When g ^ constant, the weight is modulated by the 
barycentric position (X(t)) of the population via the 
g-dependence. 

(ii) Agents with velocity i>_. An agent at position x 
and with velocity V— changes spontaneously to u+ 
with constant rate /?. This agent is allowed to ob- 
serve the environment in the observation interval O = 
[x — U, x + V] . The /3-rate is enlarged in a weighted 
proportion to the number of v+ agents present in O. 
When g ^ constant, the weight is modulated by the 
barycentric position (X(t)) of the population via the 
^-dependence. 

The following remarks can be pointed out: 

(1) Agents interact non-locally with their neighbors in 
two ways: (a) via the finite extension of the O interval 
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and (b) via the g- modulation which depends on the 
barycentcr (X (t)}. Hence, for O ^ R, one could find 
an inconstancy between the two mechanisms. This 
however is not so as the agents may simultaneously 
base their decisions on real-time information deliv- 
ered by physically different sensors like vision, sound, 
olfaction, etc. 

(2) Observe that the sum of the right- hand-sides of the 
two equations given in eqs. ^ is equal to zero. This 
is an expression of the continuity equation which guar- 
antees the conservation of the number of agents. 

(3) Using the rescaling v± n- v±/e, a H> a/e 2 and (3 M> 
(i/e 2 in eqs. (0), our two velocity model converges for 
e — > to a Burgers' type dynamics of the form: 



d x { /(M) + 7 



x+U 



^>(z,t) dz 



-v 



*(i,t) 



(6) 



which is a MF nonlinear Fokkcr-Planck equation for 
interacting diffusion processes driven by White Gaus- 
sian Noise, [THUD]. We exemplify in appendix A how 
to transfer our findings valid for the Boltzmann-type 
dynamics to this Burgers' type limiting model. 

We now explicitly give the solutions to the set of two non- 
lincarly coupled equations given in eqs. ^) for the follow- 
ing regimes: 

(A) Follow the immediate leaders : U = 0, V = e, with e 
infinitesimal small, and g = 1, 

(B) Follow all the leaders: U — 0, V — oo and g = 1, 

(C) Follow the leaders with barycentric modulation: U = 
0, V = oo and g = g(x - (X(t)}). 

Studying the regimes (A) and (B), we quantify the im- 
portance played by the size of the observation interval O. 
Comparing the regimes (B) and (C) we may appreciate 
how the influence of the g-barycentric modulation will give 
rise to flocking phase transitions. 



3 Discrete velocity dynamics 

Using a simple rescaling of the state and time variables 
(e.g. see [THj for details), we can without loss of generality 
set v + = — v— = 1 (and similarly, 7 = 1 in eq. ([6])) and 
rewrite eqs. ([2]) in canonical form: 

P ± {x,t)±d x P ± (x,t) = 

Ta(x,t)P + (x,t) ± P(x,t)P-(x,t). (7) 

Based on eq. Q, we now explore the above mentioned 
regimes. 



(A) Follow the immediate leaders 

In the myopic case where V — e is very small and U = 0, 
we may Taylor expand up to first order the quadrature 
in eqs. ^ and The resulting dynamics from eqs. Q 
reads as: 

P ± (x,t)±d x p ± (x,t) = 

±2eP+(x, t)P- (x, t) =F aP+ (x, t) ± (x, t). (8) 

This is an exactly solvable, discrete Boltzmann-type equa- 
tion discovered by Th. Ruijgrok and T. T. Wu, [25] . 
Potential applications for this type of dynamics have been 
recently considered in [TSJ [T!5] . Using the boundary condi- 
tions limui-^oo P±(x,t) = 0, the solution reads as: 

P+O, t) = -0 + d t log H(x, t) - d x \ogH{{x, t), (9) 
P-(x,t)= a-dt\ogH{x,t)-d x logH(x,t), (fO) 

where H (x, t) solves the linear Telegraphist equation: 

d u H(x, t) - d xx H{x, t) - apH(x, t) = (11) 

and whose explicit solution is recalled in the appendix B. 
The diffusive solution does not converge towards a finite 
stationary density, implying that, for any initial condition, 
the agents will ultimately be spread over the whole line 
K. This explicitly indicates that local interactions are not 
strong enough to create a cooperative motion (i.e. no flock 
is formed). 



(B) Follow all the leaders 

Instead of the myopic regime (A) given in eqs. 
now consider case (B), namely: 



), let us 



d t P ± {x,t)±d x P ± {x,t) = =p 



P-dz 



P+dz 



P + (x,t) 
P_{x,t). (12) 



(13) 



Introducing the notations 

roc 

F±(x,t)= / P ± (z,t)dz, 

J x 

we can rewrite eqs. (fl"2"|) as: 

d xt F±(x,t)±d xx F±{x,t) = ±d x [F_{x,t)F+{x,t)] 

Tad x F + (x,t)±pd x F_(x,t). (14) 

After integration with respect to x, eqs. (jT¥J) exhibit the 
same structure than eqs. §8§ for the fields F + (x,t) and 
F-(x,i). Here however, the boundary conditions are: 

lim F±(x,t) = 0, and lim F±(x,t) = p± (15) 



and the normalization imposes that (p + + pJ) = 1. With 
these boundary conditions, the solution reads as, [25j : 



F±(x,t) = <{ 1-tanh 



x — wt 



(16) 
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with 

2p ± = 1 ± w > 0, 
and where w is the velocity defined by : 



w = ^[a + p-lf -{a + /3). 



(17) 



(18) 



From cq. (TIB)) , the probability densities in eqs. (jT2J) exhibit 
the form of solitary waves: 



P±(x,t) = -d x F ± (x,t) 



P± 



8 cosh 



2 / x—wt ^ 



(19) 



The solutions in eqs. (|T^|) explicitly show that in case (B), 
the long-range imitation process generates stationary, fi- 
nite probability densities. Hence contrary to what hap- 
pens in the myopic regime (A), here flocking results from 
the long-range of the agents' mutual interactions. 

(C) Follow the leaders with barycentric 
modulation 

The above analysis of the regimes (A) and (B) suggest 
that for a certain critical interaction strength, one should 
be able to observe a phase transition from a purely diffusive 
regime with no stationary patterns to a stationary disper- 
sive flocking regime. A very simple, yet exact solution to 
this problematic is shown here, where we now introduce 
an explicit barycentric modulation function. Recent con- 
tributions using barycentric interactions are discussed in 
[H [TQ1 [T51 [21] . The class of models we consider here are of 
the following form: 

P± (x, t)±d x P ± {x,t) = +P+ {x, t) [a - F_ {x, t)] 

±P.(x,t)[/3 + ¥ + (x,t)}, (20) 



with 



W±(x,t) 



g[ Z -(X(t))}P ± (z,t)dz 



(21) 



and where (X(i)), defined in cq. (O, is the baryecnter 
position at time t. 

Let us from now on consider traveling wave solutions to 
eqs. I|20p of the form f(x — wt) for some velocity w. In such 
a traveling wave regime and for large t, we have (X(t)) = 
w t. Introducing the variable £ = (x — w t) we may, by 
a slight abuse of notations, rewrite P± in this regime as 
P±(x,t) = P±(£) and eq. §5§ now reads as: 



0= f e[p + (o+p_(o]de, 

JR 



(22) 



meaning that in the traveling wave case, the stationary dis- 
tribution, as seen from the center of mass, has zero mean. 
As the right-hand-sides in eqs. (|2H)) sum to zero, we imme- 
diately find, by simple integration with respect to £, the 
following relation: 



(l-w)P + (0 = (l + w)P-(0 + 



(23) 



where the celerity of the center of mass w will be chosen 
below, so as to be consistent with the agent distribution. 
The integration constant k will be set to zero in order 
to match natural (i.e. vanishing) boundary conditions at 
infinity. Therefore, we shall have: 



P-(0 



1 — w 
1 + w 



(24) 



Using ca. ([M| and writing eqs. (|2"0)) with the variable £, we 
obtain the MF evolution equation: 



5 ? P+(C) = 



p+(0 



1 



5 (C)P+(CK 



l 



i 



(25) 



In view of our analytical objective we now introduce the 
class of symmetric interaction modulations: 

g(x) = Acosh^a;), (26) 

with A > and 77 G R, together with the Ansatz: 

P + (0=AA(m)cosh"m (27) 

where m < is some negative constant and where Af(m) 
stands for the probability normalization, explicitly : 

v^r(M/2) 



Af{m) 



2Mr[(|m| + l)/2] 



(28) 



Introducing eqs. (j26)) and ([27)l into the integral equa 
tion (|25)l . an elementary calculation shows that eq. ([27] 
is actually a solution provided we impose: 



m + 77 = —2, 
in which case we end up with the relation: 



2AAf(m) 



1 



tanh(0 



2AM(r 



1 



w 



(29) 



(30) 



Clearly, eq. (|30p is realized only if the following equalities 
hold 



2A 



7V(r 



1 



and 



1 



1 + w 13 + 2AAf{m) '' 



(31) 



implying the self-consistency equation for the traveling 
speed: 

2aAA/~(m) 

W = 1+ m [P + 2AM(m)} (32) 

Hence, given the positive input modeling data a, (3 and A, 
together with the modulation rj of the barycentric interac- 
tion g, eq. ([2"9"]l determines the probability density decay m 
and eq. ([32]) determines the resulting solitary wave veloc- 
ity w. Note that the normalization factor in eqs. ([27l [28| 
exists only for m < 0, or equivalently for 77 > —2, which 
yields the critical power ?7 C for the range of the modulation 
decay: 

T] > r? c = -2. (33) 
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The condition on the modulation range power r\ in cq. (|33|) 
leads to conclude that for strongly localized interactions, 
r) < rjc no finite stationary probability density exists, and 
this indicates that no stable flocking does emerge from 
the agent interactions. Conversely, when 77 > rj c , a finite 
stationary probability density exists, thus showing that 
comparatively long-range interactions ultimately drive the 
population to a persistent flocking behavior. 

4 Conclusion and perspectives 

The models we have exposed in this contribution simul- 
taneously involve two types of nonlinear sources: first 
an imitation interaction mechanism of the quadratic type 
which reflects the observation of the state of the agents 
by their fellows and secondly a barycentric modulation for 
the strength of the interactions. The possibility to de- 
rive an exactly solvable model for the flocking transition 
dynamics relies partially in the fact that one nonlinearity 
can be removed. Indeed, the quadratic nonlinearity can be 
"removed" by using a logarithmic transformation, but the 
second type of nonlinearity remains as it is truly intrinsic. 
In actual agent models like birds or fishes, interactions are 
characterized by both of these nonlinearities: by imitation 
(stylized here by the quadratic nonlinearity) and by tuning 
the observation range according to the agents' local den- 
sity (effectively stylized here by a barycentric modulation). 
Another important simplification allowing the derivation 
of exactly solvable models is the mean-field limit adopted 
here. In this limit, the law of large numbers reduces the in- 
fluence of the fluctuations to insignificancy and this leads 
to deterministic evolution equations for the agent distri- 
bution. For many applications in perspective, finite-size 
population effects will manifest via the presence of fluc- 
tuations (i.e. a mesoscopic description). Noise will def- 
initely affect the dynamics and could potentially destroy 
the flocking capability of the agents. We hope that the 
basic models discussed in this paper could provide some 
analytical clues for this truly challenging issue. 

Appendix A 

The complete program performed in this paper can be re- 
peated for the Burgers' dynamics given in eq.©- Let us 
only focus here on the interactions given in case (C) for 
which the dynamics reads as: 



8t9(x,t) = -d x 



-¥(*,*) / g{z- (X{t)))y(z,t)dz 

(34) 



where: 



(X(t)) 



xV(x,t)dx. 



(35) 



Concerned with traveling wave solutions with constant 
celerity w, we introduce the notation £ = (x — wt). Ac- 
cordingly, cq. (|34p takes the form: 



*(0 



g{z)^(z)dz 



D 



a c *(0 



and w is implicitly determined by the equation: 



JW. 



(36) 



(37) 



Ea.(|36p can be reduced to the form: 

- d t log [*(£)] = j«> - J™ g (z) *(z)dz| + Ki (38) 

where k is a constant which will be taken to be zero as 
no secular probability flow exists in the stationary regime. 
Assume now a symmetric barycentric modulation of the 
form: 

g(x) = A cosh* 7 (a;), (39) 

with A > a constant and n € R. In view of ea.f|55|). we 
introduce the Ansatz: 



tf(0 =7V(m) cosh m (£), 



(40) 



where again, for m < 0, the normalization factor 7V(m) is 
given by eq. {25]). Note from eq. flU) that = #(-£) 
and therefore eq. ([57)1 is automatically satisfied. Using 
eqs. (001) and (J25J) into cq. ([351), we end up with: 



— mtanh(£) 



w-AAf(m) J cosh" +m (£)d£ 



(41) 



Using the identity J cosh(a;) 2 dx = 1 — tanh(a;), we see 
that eq.([4"Tj) is solved provided we have: 



2A 



N{m) D 



—— and w — AjV(m). 



(42) 



From eqs. (|42|) . one concludes that for 77 > —2, a traveling 
solitary wave with celerity w is created via the agents' 
interactions. Hence for 77 > —2 the flocking mechanism is 
triggered. Conversely, for r\ < —2, normalization cannot 
be achieved and this shows that too weak interactions at 
long-range preclude the formation of a flock. 

Appendix B 

The solution H(x 7 t) to the linear Telegraphist equation, 

d tt H(x, t) - d xx H(x, t) - apH (x, t) = (43) 
is of the form, [T7H25]: 

H(x, t) = - [A(x + t) + A(x - t)} 



+ ±Bx(x,t) + ^B 2 (x,t), (44) 
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where v = and where we have the following defini- 
tions: 



B 1 {x,t) = j I (y^Jf 2 - (x - z) 2 J B(z)dz, (45) 

and 

x Ij ^Vt 2 - (a; - z) 2 ) A(z)jfe, (46) 

with Iji(-) being integer-order modified Bessel functions of 
the first kind and where B(-) and A(-) are short for: 

B(x) = \ [P Q (x) - Qo(x) +a + /3] A(x), (47) 

and 

A(a:)=cxp{-±jf [P o (z) + Q o {z)-a + 0\dz}. (48) 
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